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ABSTRACT 

We provide a quantitative description and statistical interpretation of the optical 
continuum variability of quasars. The Sloan Digital Sky Survey (SDSS) has obtained 
repeated imaging in five UV-to-IR photometric bands for 33,881 spectroscopically con- 
firmed quasars. About 10,000 quasars have an average of 60 observations in each band 
obtained over a decade along Stripe 82 (S82), whereas the remaining ~25,000 have 2-3 
observations due to scan overlaps. The observed time lags span the range from a day 
to almost 10 years, and constrain quasar variability at rest-frame time lags of up to 
4 years, and at rest-frame wavelengths from lOOOA to 6OOOA. We publicly release a 
user-friendly catalog of quasars from the SDSS Data Release 7 that have been observed 
at least twice in SDSS or once in both SDSS and the Palomar Observatory Sky Sur- 
vey, and we use it to analyze the ensemble properties of quasar variability. Based on 
a damped random walk (DRW) model defined by a characteristic time scale and an 
asymptotic variability amplitude that scale with the luminosity, black hole mass, and 
rest wavelength for individual quasars calibrated in S82, we can fully explain the en- 
semble variability statistics of the non-S82 quasars such as the exponential distribution 
of large magnitude changes. All available data are consistent with the DRW model as 
a viable description of the optical continuum variability of quasars on time scales of 
~5-2000 days in the rest frame. We use these models to predict the incidence of quasar 
contamination in transient surveys such as those from PTF and LSST. 
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1. Introduction 

The optical continuum variability of quasars 
has been recognized since their first optical iden- 
tification (Matthews & Sandage 1963), and it has 
been proposed and utilized as an efficient method 
for their discovery (van den Bergh, Herbst, Prit- 
chet 1973; Hawkins 1983; Hawkins & Veron 1995; 
Ivezic et al. 2004a; Rengstorf et al. 2006). The 
observed characteristics of the variability can 
then be used to constrain the origin of their emis- 
sion (e.g., Kawaguchi et al. 1998; Trevese, Kron 
& Bunone 2001, and references therein). The 
variability of quasars has typically been quan- 
tified using a structure function (SF) analysis 
(e.g., Hughes et al. 1992; Collier & Peterson 2001; 
Bauer et al. 2009; Kozlowski et al. 2010b; Welsh, 
Wheatley, k Neil 2011), where the SF is the root- 
mean-square (rms) magnitude change (Am) as 
a function of the time lag (At) between mea- 
surements (similar to an auto-correlation func- 
tion). It is fairly well established that quasar 
variability properties depend on physical proper- 
ties such as the quasar luminosity, wavelength, 
time scale, and the presence of radio emission. 
However, despite the considerable observational 
effort invested over last few decades, many con- 
flicting claims about various correlations exist in 
the literature (see Giveon et al. 1999 for a de- 
tailed discussion). 

The traditional method for studying variabil- 
ity has been to monitor a small, select sam- 
ple of quasars over a long time baseline (e.g., 
Hawkins 2002; Giveon et al. 1999; Rengstorf et 
al. 2004). In this case, it is possible to com- 
pute the SF for each quasar, which can later be 
sample-averaged or studied individually. An al- 
ternative, utilized in more recent studies based 
on Sloan Digital Sky Survey (SDSS; York et 
al. 2000) and Palomar Observatory Sky Survey 
(POSS; Minkowski & Abell 1963) data, is to com- 
pute a single SF for all quasars in a particular 
wavelength or luminosity range. This approach, 
mandated by the fact that typically only a few 



epochs were available per object, only measures 
ensemble properties and assumes that all quasars 
selected from a narrow range of physical prop- 
erties vary in the same way. Nevertheless, this 
approach enabled studies of quasar optical vari- 
ability based on tens of thousands of objects and 
several hundred thousand photometric observa- 
tions, as well as explorations of the long-term 
variability (Vandcn Berk et al. 2004, hereafter 
VB04; Ivezic et al. 2004c, hereafter 104; Wil- 
hite et al. 2005; Mahabal et al. 2005; De Vries, 
Becker, White, Sz Loomis 2005, hereafter dV05; 
Sesar et al. 2006, hereafter Scs06). For exam- 
ple, the size and quality of the sample analyzed 
by VB04 (two-epoch photometry for 25,000 spec- 
troscopically confirmed quasars) allowed them to 
constrain how quasar variability in the rest frame 
optical/UV regime depends upon rest frame time 
lag (up to ~2 years), luminosity, rest wavelength, 
redshift, the detection of radio or X-ray emis- 
sion, and the presence of broad absorption line 
systems. By comparing SDSS and POSS mea- 
surements for ~20,000 quasars spectroscopically 
confirmed by the SDSS, Ses06 constrained the 
optical quasar variability on time scales from 10 
to 50 years (in the observer's frame). They re- 
port that there is a characteristic time scale of 
order 1 year in the quasar rest frame beyond 
which the SF flattens to a constant value. The 
SDSS has also facilitated both individual- and 
ensemble-based approaches by providing a large 
multi-epoch sample of quasars over the North- 
ern Galactic Cap and well-sampled light curves 
in the Southern Stripe 82 (S82) survey. It is re- 
assuring that the two approaches lead to simi- 
lar SFs, as discussed by dV05. A test of this 
assumption is also described in MacLeod et al. 
(2008), who show that indeed the mean behavior 
is the same. With such large samples, the en- 
semble SF(At) slopes are well-constrained, and 
the values suggest that accretion disk instabil- 
ities are the most likely mechanism causing the 
observed optical variability (VB04; Kawaguchi et 
al. 1998; see also Lyubarskii 1997). However, at- 
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tempts to constrain physical models using the en- 
semble SF are invalid as soon as one realizes that 
the ensemble SF(At) is a weighted sum of indi- 
vidual quasars with different structure functions 
(MacLeod et al. 2008). 

While studies have traditionally examined "non- 
parametric" statistical measures of variability 
such as the SF, a major challenge has been 
to describe the variability of individual quasars 
in a compact way. Recently, the introduction 
of a damped random walk (DRW) model has 
provided a way to mathematically characterize 
quasar light curves in terms of a characteristic 
time scale (r) and an amplitude (SFoo) which 
are then correlated with the physical properties 
such as luminosity and black hole mass. Kelly et 
al. (2009, hereafter KBS09) modeled a sample of 
100 quasar light curves as a DRW and suggested 
that thermal fluctuations driven by an underly- 
ing stochastic process such as a turbulent mag- 
netic field may be the dominant cause for the op- 
tical flux fluctuations. Kozlowski et al. (2010a; 
hereafter KozllO) applied the DRW model to 
the well-sampled Optical Gravitational Lensing 
Experiment (OGLE) light curves (Udalski et 
al. 1997; Udalski et al. 2008) of mid-infrared- 
selected quasars behind the Magellanic Clouds 
from Kozlowski & Kochanek (2009). Their anal- 
ysis shows that the DRW model is robust enough 
to efficiently select quasars from other variable 
sources, despite the large surface density of fore- 
ground Magellanic Cloud stars (see also Butler 
& Bloom 2011; MacLeod et al. 2011; Kozlowski 
et al. 2011a). 

In MacLeod et al. (2010, hereafter MaclO), 
we applied the DRW model to the light curves 
of ~ 10,000 quasars in S82 and found a correla- 
tion between SF^o and black hole mass which is 
independent of the anti-correlations with lumi- 
nosity and wavelength (see also Ai et al. 2010; 
Meusinger et al. 2011). We also found that 
T increases with increasing wavelength, remains 
nearly constant with redshift and luminosity, and 
increases with increasing black hole mass (see 



also KBS09; KozllO). In Kelly et al. (2011), it 
was shown that a similar stochastic model but 
with multiple time scales for a single object can 
accurately reproduce the X-ray variability of ac- 
tive galactic nuclei (AGN) and microquasars. An 
inhomogeneous accretion disk model, where the 
temperature fluctuations throughout the disk are 
driven by a DRW process, can explain the disk 
sizes derived from microlensing light curves (see 
Morgan et al. 2010) while matching the observed 
level of optical variability, and predicts SEDs 
which are in better agreement with observations 
than standard thin disk models (Dexter &; Agol 
2011). 

One defining feature of the DRAV model for a 
single quasar is that it predicts a Gaussian distri- 
bution of magnitude differences Am for a given 
At. On the other hand, the observed Am distri- 
bution in the optical for an ensemble of quasars 
observed at two times separated by At deviates 
strongly from a Gaussian but is well fit by an ex- 
ponential distribution (104). This conflict repre- 
sents an important puzzle for understanding the 
DRW model and its applicability to quasar light 
curves. Also, the high likelihood of extreme val- 
ues of Am has important implications for the 
interpretation of observations of transients. For 
example, Vanden Berk et al. (2002) reported the 
detection of an orphan gamma-ray burst after- 
glow based on the 2.5 magnitude decrease in op- 
tical flux. Such a large flux change was inconsis- 
tent with a quasar based on a Gaussian model 
for their variability, but the source was never- 
theless confirmed to be a highly variable quasar 
(Gal- Yam et al. 2002). An accurate statistical 
description of the two-epoch photometry for en- 
sembles of quasars will be important for transient 
detection in large surveys, where quasars repre- 
sent a major contaminant. 

Our goal here is to produce a unified view of 
ensemble and individual optical variability in the 
context of the DRW model. In this study, we 

show that the differences in shape between the 
ensemble SF and the DRW SF are well-explained 
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by averaging over the properties of individual 
quasars. We also show that the exponential dis- 
tributions of magnitude changes for ensembles 
of quasars at fixed time lag are naturally con- 
structed by summing the intrinsically Gaussian 
distributions of magnitude changes produced by 
individual quasars. There are several residual is- 
sues which we discuss as part of the compari- 
son. An overview of the SDSS and POSS data 
used in this study is presented in Section 2. In 
Section 3, we describe the observed ensemble 
quasar variability in terms of the SF as a func- 
tion of wavelength and time lag in the observer's 
frame. We then convert to rest-frame quantities 
and compare the data to a model ensemble SF 
based on our previous DRW analysis of S82 light 
curves. In Section 4, we combine the constraints 
on short-term quasar variability based on SDSS 
data with the constraints on long-term variabil- 
ity derived from matching the SDSS and POSS 
catalogs. In Section 5, we discuss the implica- 
tions our results have on transient identification, 
with a focus on future time-domain surveys such 
as the Palomar Transient Factory and the Large 
Synoptic Survey Telescope. Our results are dis- 
cussed and summarized in Section 6. 

2. Data Overview 

In this section, we briefly summarize the rel- 
evant SDSS and POSS data. We focus our de- 
scription on the quasars with multiple observa- 
tions (similar to MacLeod et al. 2010; 2011). 

2.1. The Basic Characteristics of the 
SDSS Imaging Survey 

The SDSS provides homogeneous and deep 
(r < 22.5) photometry in five passbands {u, g, 
r, i, and z, Fukugita ct al. 1996; Gunn et al. 
1998; Smith et al. 2002; Ivezic et al. 2004b) 
accurate to 0.02 mag, of up to 12,000 deg^ in 
the Seventh Data Release (DR7, Abazajian et 
al. 2009). The DR7 sky coverage results in 
photometric measurements for about 357 million 



unique objects. Astrometric positions are ac- 
curate to better than 0.1 arcsec per coordinate 
(rms) for sources brighter than 20.5 mag (Pier 
et al. 2003), and the morphological information 
from the images allows robust star-galaxy sep- 
aration to ~21.5 mag (Lupton et al. 2001). A 
compendium of technical details about the SDSS 
can be found in Stoughton et al. (2002). 

The SDSS offers an unprecedented photomet- 
ric accuracy for such a large scale optical sur- 
vey. Not only arc the photometric errors gen- 
erally small, but they are accurately determined 
by the photometric pipeline {photo, Lupton et 
al. 2001) and can be reliably used to estimate 
the statistical significance of measured magni- 
tude differences. This ability is of paramount 
importance for a robust statistical study of vari- 
able objects. This error behavior is illustrated in 
Figure 2 in Ivezic et al. (2003). Throughout our 
analysis, we assume the SDSS photometric errors 
to be CTphot = 0.018 mag in g, r, and i bands, and 
0.04 mag in u and z (sec also Figure 1 in Sesar et 
al. 2007). The photometric errors are assumed to 
be independent of magnitude due to the bright 
sample limit oi i < 19.1 (see next section), and 
for simplicity when accounting for photometric 
errors in the SF calculations. 

2.2. The SDSS Multi-epoch Data 

The SDSS imaging data are obtained by imag- 
ing the sky in six parallel scanlines, each 13.5 ar- 
cmin wide (a "strip" in SDSS terminology). The 
six scanlines from two adjacent scans are then 
interleaved to make a filled "stripe". Because 
of the scan overlaps, and because the scans con- 
verge near the survey poles, ^^40% of the sky in 
the northern survey is essentially surveyed twice. 
In addition, 290 deg^ of the southern survey area 
lies along Stripe 82 (S82) and has been observed 
about 60 times on average to search for vari- 
able objects and, by stacking the frames, to go 
deeper (Frieman et al. 2008; Annis et al. 2011). 
This valuable subsample contains well-sampled 
light curves for 9,258 spectroscopically confirmed 
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quasars whose variability properties are analyzed 
in MaclO, and it can be used to verify some of 
the results inferred from the analysis of two re- 
peated observations. Overall, the SDSS has ob- 
tained multi-epoch data for ~4,000 deg^ of sky, 
with time scales ranging from 2 hours to over 9 
years, and with a wide range of Galactic latitudes 
extending all the way to the Galactic plane. 

We define a quasar as any object listed in 
the SDSS catalog of spectroscopically confirmed 
quasars (the "DR7 Quasar Catalog", Schneider 
et al. 2010). Its most recent fifth edition lists 
the SDSS DR7 BEST photometry for 105,783 
quasars from 9,380 deg^. For a description of 
the spectroscopic target selection for quasars, see 
Richards et al. (2002). We note that the qual- 
ity of photometry is good for the objects in the 
DR7 Quasar Catalog. In total, there are 33,881 
spectroscopically confirmed quasars with at least 
two observations. We provide a catalog of all 
the SDSS repeated imaging of quasars online^^ 
(see the Appendix for a detailed catalog descrip- 
tion). The sky distribution of repeatedly imaged 
quasars is shown in Figure 1. 

When we compare pairs of observations, we 
use the secondary imaging observation with the 
largest time lag from the primary imaging ob- 
servation listed in the DR7 Quasar Catalog, for 
a total of two observations per quasar. We de- 
fine Am = 772-2 — n^-i, where m2 is the magni- 
tude observed at a later epoch than mi. We 
omit the S82 quasars since we do not wish to 
test our model on the same sample of objects 
from which the DRW model parameters were 
derived. The DR7 Quasar Catalog is spectro- 
scopically complete to i < 19.1 in the quasar re- 
gion of color space. Therefore, unless otherwise 
noted, we require that the primary observation 
for the remaining 24,627 quasars has i < 19.1 
in order to be consistent with the flux limit of 
the quasar sample, for a total of 14,939 quasars. 
Given the SDSS limit of ~22 mag, we further re- 
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Fig. 1. — Distribution of 33,881 DR7 quasars with at 
least 2 observations on the sky in equatorial coordinates. 
The number of observations increases towards the survey 
poles where the "stripes" overlap. The inset shows a closer 
view of the SDSS Stripe 82, which has ~60 observations 
per object (the density of quasars is non-uniform in S82 
due to the increasing Galactic contribution toward nega- 
tive R. A.). 



strict the magnitude differences in each band to 
I Am I < 3 mag in order to reduce the amount of 
contamination by poor photometry, for a total of 
~80,000 pairs of SDSS measurements (summed 
over all bands). Finally, we adopt the redshifts 
and absolute magnitudes listed in the Schneider 
et al. (2010) spectroscopic quasar catalog, and 
the black hole masses as measured from emission 
lines by Shen et al. (2011). All magnitudes are 
corrected for Galactic extinction using Schlegel 
et al. (1998). 

2.3. SDSS POSS Long-term Measurements 

While SDSS S82 obtained measurements with 
time differences of up to 10 years, longer time 
scales of up to 50 years can be probed by com- 
paring the SDSS and POSS catalogs. Ses06 have 
addressed the problem of large systematic er- 
rors in POSS photometry by recalibrating several 
pubhcly available POSS catalogs (USNO-A2.0, 
USNO-Bl.O, DPOSS and GSC2.2). A piecewise 
recalibration of the POSS data in 100 arcmin^ 
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patches (one SDSS field) generally resulted in 
an improvement of photometric accuracy (rms) 
by nearly a factor of two compared to the origi- 
nal data (POSS-I magnitudes can be improved to 
~0.15 mag accuracy, and POSS-II magnitudes to 
~0.10 mag accuracy). In addition to the smaller 
core width of the error distribution, the tails of 
the distribution become much steeper after the 
recalibration. These improvements are mostly 
due to the very dense grid of calibration stars 
provided by the SDSS, which rectifies the intrin- 
sic inhomogeneities of Schmidt plates. 

The much longer time lags between the Palo- 
mar Observatory Sky Surveys (POSS-I and POSS- 
II) and the SDSS (up to ~50 years) than spanned 
by the available SDSS data make it easier to de- 
tect deviations of the SF from a simple power- 
law. dV05 and Ses06 compared the SDSS and 
POSS data for over 10,000 quasars from the 
SDSS Data Release 2 in order to constrain the 
long-term quasar variability. As discussed by 
Ses06, the dV05 and Ses06 measurements of the 
SF agree within a la uncertainty of the Ses06 
measurements. Here, we use results from both 
the dV05 and Ses06 studies. We also make use 
of data from the Digitized Palomar Observatory 
Sky Survey (DPOSS; Djorgovski et al. 1998) that 
overlap 8,000 deg^ of sky from the SDSS Data 
Release 5. These DPOSS data were recalibrated 
following the procedure outlined by Ses06. There 
are 81,189 SDSS DR7 quasars with DPOSS ob- 
servations in the G, R, and / bands, and we 
provide their 2-epoch photometry at the same 
website (see footnote 12). For our analysis in 
Section 4, we use the primary SDSS observations 
as listed in the DR7 quasar catalog when com- 
paring SDSS-DPOSS magnitudes. We calculate 
synthetic POSS GRI magnitudes from the SDSS 
photometry, require i < 19.1, and impose var- 
ious quality cuts following Ses06. In total, we 
have 56,732 SDSS-DPOSS Am measurements. 



3. Characteristics of Observed Variabil- 
ity 

For our SF analysis, we examine the behavior 
in both the observer's and the quasar rest frames 
because the former is informative for data inter- 
pretation (e.g., transients in large surveys) and 
the latter constrains quasar physics. The SF be- 
havior in the observer's frame is discussed in Sec- 
tion 3.1, and we convert to rest-frame quantities 
in Section 3.2. We handle various trends with 
observed properties by taking narrow ranges of 
the relevant quantities, in particular, the time 
lag and wavelength in each frame. The luminos- 
ity, mass, and redshift information are utilized 
in Section 3.2, where we test whether the scal- 
ings with these physical parameters derived from 
individual light curves in S82 can reproduce the 
ensemble variability of the two-epoch sample. 

3.1. Quasar Variability in the Observer's 
Frame 

Assuming that the observed variability reflects 
the physics in the accretion disk, the behavior 
in the observer's reference frame will be a con- 
volution of rest-frame variability over redshift, 
luminosity, and other parameters. While this 
convolution will obscure important physical scal- 
ings, quasar variability in the observer's reference 
frame is still of major interest when interpret- 
ing survey data, in particular when distinguish- 
ing between quasars and other variable sources 
such as transients, or in using variability to se- 
lect quasars. Therefore, we start by considering 
quasar variability in the observer's frame. 

3.1.1. Am Distribution 

Figure 2 shows the cumulative distribution of 
Am in the SDSS u, r and z bands for four dif- 
ferent At ranges. The thin lines show the pre- 
dicted distributions based on Gaussian (dotted) 
and exponential (dashed) analytic functions with 
the same rms as the data. The photometric er- 
rors are taken into account by adding a Gaussian 
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Fig. 2. — The observed cumulative distributions of urz 
magnitude differences as a function of time lag in the ob- 
server's frame. We define Am = m2 — mi, where m2 is the 
magnitude observed at a later epoch than mi. The data 
(thick lines with Poissonian error bars) are only shown 
for sample sizes exceeding 10 data points (the probabili- 
ties become increasingly unreliable when below ~0.003). 
Going from the narrowest distribution to the widest, the 
time lags span 1-30 days, 50-150 days, 200-400 days, and 
1400-3000 days. The thin lines show the predicted dis- 
tributions based on Gaussian (dotted) and exponential 
(dashed) analytic functions with the same rms as the data 
for each time lag range, where the photometric errors are 
taken into account by adding a Gaussian component of 
width defined in Section 2.1. 

component of width aphot = 0.018 mag in r and 
0.04 mag in u and z. The exponential curves pre- 
dict a much higher probabihty of large magnitude 
changes than the Gaussian curves, and the data 
fohow this prediction (we only show the data for 
bins with more than 10 points). The data points 
become increasingly unreliable in the tails of each 
distribution due to small sample sizes, and this 
fact may cause the large discrepancy with the 
exponential curve for small At in the r and z 
bands. 

Figure 3 shows the differential Am distribu- 
tions for each SDSS band and for two slices of 
observed-frame time lag. The distribution width 
(i.e., the SF value) increases with time lag and 
decreases with wavelength. The wings of the dis- 
tributions are also closer to exponential (more 
accurately, a double-exponential or Laplace dis- 
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Fig. 3. — The symbols show the distribution of mea- 
sured magnitude differences. Am, in the five SDSS bands, 
and for two narrow ranges of time lag. At, as indicated at 
the top. n indicates the number of points in a bin divided 
by the bin width, and Ntot is the total number of points 
used for each histogram. The distribution width (a) calcu- 
lated using Eq. 1, shown in each panel, is increasing with 
time and decreasing with wavelength. The dot-dashed 
lines show Gaussian distributions with a root-mean-square 
(rms) equal to a, and the dashed lines show exponential 
distributions, exp(— Am/Ac), with Ac = rms/V2. Note 
that an exponential distribution provides a better fit to 
the wings of the observed distributions. Each panel also 
displays the number of measurements outside the ±3a 
range (A'^ > 3a), their corresponding fraction (/) ex- 
pressed in percent (the expected value of this fraction for 
a Gaussian distribution is 0.3%, and for an exponential 
distribution 1.4%), and the ratio R — rms/o". 



tribution) than to Gaussian. Figure 4 shows sim- 
ilar distributions for the changes in color for two 
slices in time lag. The color changes are smaller 
than for the individual magnitudes at the same 
time separation. In the bottom six panels, the 
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color changes are plotted against each other as 
well as against the r-band differences. Although 
the scatter around these relationships is quite 
large, quasars tend to get bluer as they brighten 
on average (bottom four panels), in accordance 
with previous results (e.g., Giveon et al. 1999). 
An analogous correlation was found recently by 
Schmidt et al. (2011), who analyzed individual 
light curves. Their study showed the slope Sgr 
between the r-band and 5-band variations, as in 
r — (r) = {sgr + l){g — (g)) + b, is on average 
around —0.2, indicating that quasars get bluer 
as they brighten. This value corresponds to a 
slope of 4 in the bottom-right panel of Figure 4 
(equating m — (m) to Am), shown by the solid 
line, which approximately follows the observed, 
two-epoch distribution. 

3.1.2. Structure Function 

Next, we quantify the quasar variability using 
a SF analysis. Although several definitions can 
be found in the literature, the SF essentially mea- 
sures the rms magnitude difference function 
of time lag between magnitude measurements. 
We compute the structure function as 

SF = 0.74(IQR)/ViV^, (1) 

where IQR is the 25% - 75% interquartile range 
of the Am distribution, and N is the number of 
Am values. This approach is insensitive to out- 
liers in the data which may result from poor data 
quality and is especially effective at short time 
lags where the SF is small. This value is equiv- 
alent to the rms if the distribution is Gaussian. 
When stated, the photometric errors are taken 
into account by subtracting V^(Jphot in quadra- 
ture from the SF. Two other definitions of the SF 
are found in the literature. Following Bauer et 
al. (2009), we refer to them as SF^^) and SF^^); 

SF(^)(Ai) = ^{Am^) (2) 



SF(^)(At) = ^-(|Am|)l (3) 



Left: 300cl < At < 400cl Right: 1400d < At < 3000d 
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Fig. 4. — The top three rows are analogous to Fig. 3, 
except that we show the changes in the u—g, g—i, and i—z 
colors for two time lags in the observer's frame (left: 300- 
400 days, right: 1400-1600 days). The fourth row shows 
the change of u — g color as a function of the change of 
r — i color, and the bottom two rows shows the change of 
the r band magnitude as a function of the change of <? — i 
OT g — r color. The contours show regions containing 25%, 
50%, 75% and 90% of the data points. The solid line in 
the bottom-right panel shows the slope found by Schmidt 
et al. (2011) using individual quasar light curves. 

For a Gaussian distribution, the ratio of the 
form adopted here to the other two forms is 
SF/SF(^) = SF/SF(-^) = 1, while for an ex- 
ponential (Laplace) distribution, SF/SF^"^^ = 
1.03/^/2 and SF/SF^-^) = 0.82. 

Figure 5 compares the SF for the combined 
SDSS and POSS data over time scales ranging 
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Fig. 5. — The structure function for quasar variability measured in the g band and in the observer's frame (corrected 
for errors using SF — ^ySF^^^^ — 2(7^^^^). The red data points show the SDSS measurements, and the data points at 
At > 3000 days are inferred from comparing SDSS data and POSS I/II data, taken from Ses06 (open squares) and de Vries, 
Becker, and White (2003, blue triangles, which have similar errors as those from Ses06). The short-dashed line shows the 
best power-law fit to the SDSS measurements alone, SF oc At^, and the long-dashed line shows a simultaneous fit to all 
data points, SF — Saa[i — exp(— Af/rsp)]^''^. The best-fit parameters for these two fits are listed in the upper-left corner. 
The dotted line shows the prediction of the DRW model trained on S82, as described in Section 3.2.1. 



from 5 days to 50 years in the observer's frame. 
The SF values shown up to At ~ 2000 days 
are computed using Eq. 1, with the photomet- 
ric errors subtracted. For variability in the g 
band with time lags of 10 and 50 years, we 
adopt 0.35 and 0.43 mag to represent the dV05 
results, and 0.21 and 0.40 mag for the Ses06 
results, respectively. The errors in these esti- 
mates are ~0. 05-0. 10 mag. The short-dashed 
line shows the best power-law fit to the SDSS 
measurements alone (SF oc At*^'^^). It over- 
estimates the variability level for time scales 
longer than several years. The long-dashed line 
shows the best-fit asymptotic function to all data 
points given by Soo [1— exp(— At/rsp)]^^^) a com- 
mon parametrization of the SF (e.g., Hook et 
al. 1994), with Soo and tsf as free parame- 
ters. The limiting value of the best-fit over- 
all structure function, 5oo = 0.28 ± 0.01 mag, 
is probably systematically uncertain at a level 



of 0.02-0.03 mag due to the discrepancies be- 
tween the SDSS-POSS-I/II measurements (see 
Ses06). The measured characteristic time scale, 
Tsf ~ 1400 days, corresponds to ~700 days in the 
rest frame and is uncertain by ~10-20%. While 
this analysis and the earlier ones by dV05 and 
Ses06 provide strong evidence that the SF lev- 
els out on long time scales, we cannot rule out a 
continuing but slower rise. 

3.2. Quasar Variability in the Rest Frame 

The redshift distribution of the quasars en- 
ables one to map a discrete distribution of wave- 
lengths and time differences At in the observer's 
frame to a smoother distribution in the quasar 
rest frame. Figure 6 shows the distribution of 
our 79,787 Am measurements for i < 19.1 ob- 
jects in the rest-frame time difference {AtRp) 
and wavelength (Xrp) plane. Most (54%) ob- 
servations have Atjip < 50 days. The discrete 
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distribution^^ of Ao and Ato in the observer's 
frame is spread along the hnes Atjip = kXjip, 
with k = Ato/Xo, and according to the sam- 
ple redshift distribution [Atjip = Ato/(l + z), 
A_RF = Ao/ (l+z)]. Here, we consider the shape of 
the Am distribution as a function of these quan- 
tities, the form of SF {At jip), and the dependence 
of the SF on all physical parameters simultane- 
ously, including the i-band absolute magnitude 
Mi (K-corrected to the rest frame) and the black 
hole mass Mbh- 

3.2.1. Application of a DRW model 

We develop a model SF which reproduces the 
observed ensemble variability using information 
derived from individual quasar light curves in 
S82. Our model is motivated by the success of a 
DRW in describing quasar light curves (KBS09; 
KozllO; MaclO), as well as the detection of a 
turnover in the ensemble SF on long time scales, 
suggesting a characteristic time scale for variabil- 
ity (e.g., Ses06; Welsh, Wheatley, & Neil 2011). 
We assume that the turnover is a consequence 
of the characteristic time scale distribution mea- 
sured among individual quasars, and that the 
turnover in the ensemble SF corresponds to the 
average time scale, tsf- Suppose each quasar is 
described by its own DRW parameters, r and 
SFoQ. Then, the ensemble SF is a weighted con- 
tribution of all the individual structure functions 
over their distribution in these parameters, 

SF(AO = / drdSFoo ^^^gp SF(At|r,SFoo)gso, 

(4) 

where SF(At|r, SF^) 

qso is the structure function 
at time At for a quasar with DRW variability 
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^^Note that a discretely-distributed observed At can lead to 
artificial correlations between the SF and I^Irf ■ For fixed 
Ato, longer ISIrf correspond to both lower redshift and, 
because of magnitude limits, lower luminosities. Since 
there is also an anti-correlation between variability am- 
plitude and luminosity, this leads to a SF that increases 
toward higher IStRF- Therefore, caution must be taken 
when binning in AtRp using very sparse Ato so that one 
does not mistake wiggles in the ensemble SF(At_RF) for 
multiple intrinsic time scales, for example. 



Fig. 6. — (Top) The distribution of SDSS magnitude dif- 
ference measurements in the plane spanned by rest-frame 
wavelength, Xrf and time lag, AtRp. The discrete dis- 
tribution of Ao and Ato in the observer's frame is spread 
along the lines AtRp ~ kXRp, with k = Ato/Xo, and 
according to the sample redshift distribution [AtRp = 
Ato/ {1 + z), Xrf ~ Ao/ (1 + z)]. The bottom panel shows 
the number of measurements in each bin on a log scale 
according to the legend at top. The white pixels contain 
fewer than 5 data points. 
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parameters r and SFqo, given by 

SF(At|r, SFoo),so = SFoo(l - e-l^*!/-)^^. (5) 

In particular, SF(At)gso is the expected stan- 
dard deviation of magnitude differences Am for 
a given quasar at a time lag At. SF{At)qso is 
related to the standard deviation in magnitudes 
(fjm) at a given At by SF(At)gso = \/2(Jm, where 
the factor of ^/2 results from subtracting two 
magnitudes. 

To build a model for the ensemble variability, 
we follow these steps for each quasar in the two- 
epoch sample: 

1. Predict r and SFqo based on the quasar's 
physical parameters. 

2. Include the intrinsic scatter in r and SFqo 
for quasars with similar physical parame- 
ters. 

3. Estimate the SF value at the measured 
time lag At up using Eq. 5. 

4. Draw one model Am value from a Gaussian 
distribution with a width set by SF(Ati^i?)gs 
adding photometric noise if necessary. 

We generally average over 1000 of these Monte 
Carlo models for the distributions expected from 
the DRW model. 

In the first step, r and SFqo are estimated us- 
ing the scalings of the DRW parameters with the 
quasar's physical parameters found in MaclO: 

log/ = ^ + 51og(^)+C7(M. + 23) 

where A = -0.51, B = -0.479, C = 0.131, and 
D = 0.18 for / = SFoo (in mag), and A = 2.4, 
B = 0.17, C = 0.03, and D = 0.21 for / = r 
(in days). Note that the dependence on redshift 
was found to be negligible in both cases {E = 
0). By applying this model to the quasars in 



the large two-epoch sample, we are also testing 
the accuracy of the scalings derived for the S82 
quasars. 

In the second step, we account for the intrin- 
sic scatter in r and SFoo for quasars with simi- 
lar physical parameters. When measured for in- 
dividual quasar light curves, r and SFqo show 
scatter about their mean trends (Eq. 6). The 
magnitudes of these residuals are too large to 
be fully attributed to measurement uncertainties 
(MaclO; Bauer et al. 2011). From simulations 
(see Section 4.2 of MaclO), we estimated that the 
latter statistical uncertainties account for only 
70%, 60%, and 13% of the scatter in r, SFoo, and 
(T = SFoo/ y/r, respectively. The fitting errors are 
much smaller for a because it is well-constrained 
even for light curve lengths shorter than r. If 
we define K = r-v/SF variable orthogonal 

to a (in log space), its uncertainty due to fitting 
errors contributes ~82% of its scatter. Figure 7 
compares the scatters in observed a (top panel) 
and K (bottom panel) for S82 quasars to those 
for Monte Carlo models of the light curves. 

Figure 8 shows the distributions of a and K 
measured for the S82 sample compared to their 
expected values from Eq. 6. We find that the dif- 
ferences between the observed log a (JogK) and 
the estimates from Eq. 6 peak near zero with 
an rms of 0.16 (0.62) dex. After taking the fit- 
ting errors into account, the rms is reduced to 
0.149 (0.26) dex, which we take as the intrinsic 
stochasticity. Therefore, in the second step, we 
take the r and SFoo estimated from Eq. 6, com- 
pute a and K, add a random Gaussian deviate 
of width 0.149 (0.26) dex to each logo- (logi^) 
value, and then convert back to r and SFoo. This 
process should provide a reasonable model for the 
intrinsic scatter. 

In the third step, the SF for a particular 
quasar at a given At^p is estimated using Eq. 5. 
In the fourth step, a model Am value is drawn 
from a Gaussian distribution with a standard de- 
viation of SF{AtRF)gso, the Gaussian distribu- 
tion expected from the DRW model. To account 
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Fig. 7. — Determination of fitting errors. The distribu- 
tions of DRW parameters a - 



and K = ry/WZ. 
for Monte Carlo models of S82 quasar light curves are 
shown in the top and bottom panels, respectively. The 
solid, black curve in each panel shows the distribution of 
observed g (K) values, which are used as inputs to gen- 
erate the light curves, normalized by the median value of 
0.16 mag yr"^''^ (200 mag^^^ days). The red, dashed line 
shows the ratio between the best-fit parameters for the 
simulated light curves and their input values. The red 
histogram is narrower than the black histogram since the 
resulting best fit should be similar to the input value that 
generated the light curve. The contribution of fitting er- 
rors to the overall scatter of observed values is estimated 
to be the ratio of the variances (cr^) of the red and black 
histograms. 



for photometric errors, we further add a random 
Gaussian deviate of width \/2apfiot to each model 
Am value (the factor results from adding the 
photometric errors in quadrature). This proce- 
dure results in one model Am value per quasar 
based on the expectations from the DRW model. 
We then repeat this 1000 times for each quasar in 
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Fig. 8. — Distributions of ct = SFoo/^ and K = 
T\/SFoo for S82 quasar light curves with r < 19, nor- 
malized by the expected values from Eq. 6. Each dot- 
dashed curve shows a Gaussian with unit area, with a 
mean (fi) and rms (a) as listed in the legends. After re- 
moving the contribution from fitting errors (F — 13% 
and 82% for a and K, respectively), the intrinsic dis- 
tribution widths are ai„t = 0.149 and 0.26 dex, where 

CTznt = (f^' - F<jV100)^/^ 



the sample. The resulting model for the ensem- 
ble SF is then the rms width of this model Am 
distribution otherwise calculated in the same way 
as for the data (using the interquartile range). 

3.2.2. Explaining the Exponential Tails of the 
Am Distribution 

In Figure 9, the observed magnitude differ- 
ence distribution is shown for a narrow slice in 
AtRF and Xrf- The observed distribution is 
very similar to the model distribution (solid line) , 
in which each quasar is assigned a model Am 
value drawn from a particular SF{AtRF)qso, as 
described above. The exponential distribution 
for large |Am| results from a superposition of 
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Fig. 9. — The red data points with errors bars show 
the observed magnitude difference distribution for 365 < 
AtRF < 730 days and 2000 A < Xrf < 3000 A. The sohd 
hne shows the expected distribution based on the DRW 
modeL The dashed Gaussian curves show the expected 
distributions for 9 individual quasars with different values 
of SF{AtRF)qso (and thus r and SFoo; see Eq. 5). Each 
curve has been convolved with a Gaussian noise compo- 
nent of width (Jphot to account for photometric errors, and 
shifted by the median of the observed distribution. 

many Gaussians (dashed lines) corresponding to 
different values of SF{Atjip)qso- The range in 
SF {At jip) qso is caused by differing r and SF^o 
values, which can be attributed to a range in 
quasar luminosity and black hole mass plus in- 
trinsic scatter. 

Figure 10 shows the distribution of Am in 
four bins of rest wavelength and three bins of 
rest-frame time lag. The distributions remain 
exponential at large \Am\ for all 12 combina- 
tions of AtfiF and Xrf, as illustrated by the 
dashed lines. The model distributions, which 
carry information about each individual quasar's 
expected DRW parameters (and thus incorporate 
the Mi and Mbh information) , are shown as solid 
curves. The observed and predicted SF {At rf) qso 
values (cj) are listed to the right and left of each 
histogram, respectively. The model and data dis- 
tributions agree well, showing that the exponen- 



tial distributions seen in the statistics of ensem- 
bles of quasars naturally result from averaging 
over quasars that are individually well-described 
by a Gaussian DRW process. There are two sys- 
tematic discrepancies, however. First, while the 
exponential tails can be reproduced at large rest- 
frame time lags (500 < At^p < 1500 days), the 
value of the observed SF is systematically higher. 
Since there are no known intrinsic differences in 
the physical properties of the S82 quasars and 
the two-epoch sample studied here, the discrep- 
ancy is likely due to biases in estimates of the 
DRW variability parameters for the S82 quasars. 
Secondly, at small time lags (50 < Atjip < 150 
days), the model over-predicts the data rms by 
10%, while the statistical uncertainty is 1% as- 
suming a perfect Gaussian distribution of Am 
and no other systematic errors. This discrep- 
ancy may be due to some systematic effect that 
was not accounted for when computing the error 
bars (or when correcting forthe SDSS photomet- 
ric errors). Alternatively, the discrepancy could 
result if the DRW model is inaccurate on these 
time scales (1-200 days). These discrepancies are 
discussed further in Section 4.1. 

3.2.3. The SF as a Function of Time Lag and 
Wavelength 

Now that we understand and can reproduce 
the shape of the distribution of magnitude changes 
Am between two times, we can simply consider 
the SF (i.e., the width of this Am distribution) 
as a function of physical quantities such as AtRF 
and Xrf and test the model prediction. Fig- 
ure 11 shows the SF as a function of AtRp. 
Here, the measured SF is not corrected for the 
SDSS photometric accuracy. Again, it is ap- 
parent that the model predicts a systematically 
lower SF for time lags less than 200 days. How- 
ever, the stochastic model predicts a systemati- 
cally higher SF for the longest time lags. This 
suggests the model SF may be biased low at long 
time lags, and the bias seems to be most promi- 
nent at shorter wavelengths (Figure 12). Note 
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Fig. 10. — As in Fig. 3, except that the symbols show the distribution of measured magnitude differences, Am, for 
subsamples selected using rest-frame quantities, as marked in the panels. Each row corresponds to the same wavelength 
range, and each column to the same AtRp range. The solid lines show the predicted DRW model distributions (see 
Section 3.2.1) with widths a listed to the right of the histograms. 
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Fig. 11. — The red symbols show the structure func- 
tion for quasar variabihty as a function of rest-frame 
time lag, for four ranges of rest-frame wavelength, as 
marked. The short-dashed line in each panel is the best- 
fit power law (SF oc At''), with /3 marked in each panel. 
The long-dashed lines are the best-fit exponential curves, 
SF = Soo[l — exp(— At/rsp)]"'^''^, with the parameters 
listed in each panel. Both curves include a photomet- 
ric noise component equal to the minimum measured SF 
in each panel (marked by the horizontal dot-dashed line). 
The dotted lines with cross symbols show the DRW model 
prediction (for 100 realizations, see Section 3.2.1). The 
green dotted lines with circles show the prediction of the 
DRW model trained on S82 data with all r and SFoo val- 
ues multiplied by 2 and \/2, respectively. 



that the measurements at long rest-frame time 
lags and short rest wavelengths are dominated 
by the u band, so there may be biases simply 
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Fig. 12. — The red symbols show the structure func- 
tion for quasar variability as a function of rest-frame 
wavelength, for five ranges of rest-frame time lag, as 
marked. The dashed lines are the best-fit power laws 
(SF oc A") with the power- law index listed in the top- 
left corners. The dash-dotted lines show the functional 
dependence adopted by 104 (a = —0.3). Both curves in- 
clude a photometric noise component corresponding to 
Cphot ~ 0.02 mag. The dotted lines with cross symbols 
show the DRW model prediction (for 100 realizations, see 
Section 3.2.1). The green dotted lines with circles show 
the prediction of the DRW model trained on S82 data with 
all T and SFoo values multiplied by 2 and \/2, respectively. 



from the fact that the 14% of u and z light 
curves which were dominated by noise were omit- 
ted from the analysis in MaclO. We investigate 
such biases further in Section 4.1. 

In Figure 13, we show the observed ensemble 
SF without correction for the SDSS photomet- 
ric errors. We also show the signal-to-noise ratio 
for each bin (top-right panel). In the bottom- left 
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Fig. 13. — The top-left panel displays the structure 
function for quasar variability measured by the SDSS as 
a function of rest-frame wavelength and time lag. SF is 
computed as 0.74(IQR)/VA^ - 1, where IQR is the 25% 
- 75% interquartile range. The figure is based on 79,787 
SDSS measurements, and shows the structure function on 
a linear scale from 0.03 to 0.4 mag, as indicated in the 
top bar. The top-right panel shows the SF signal-to-noise, 
where the error asF is computed as 1.15SF/y/ N — 1. The 
bottom-left panel shows the model structure function (see 
text) , and the bottom-right panel shows the (data-model) 
residuals. The red line shows the SDSS composite quasar 
spectrum in arbitrary flux (Fx) units from Vanden Berk 
et al. (2001). The weak vertical feature in the residual 
map at ~2800A is coincident with the Mgll line visible in 
the composite spectrum. The white pixels in each panel 
contain fewer than 5 data points even when including all 
adjacent bins (and thus are not used in our analysis). 



panel, we show the model expectation including 
the estimated photometric errors. The bottom- 
right panel shows the (data - model) residuals. 
Overall, the results are excellent, with a median 
difference of only 0.02 mag and a median absolute 
difference of 0.03 mag. We see a slight decrement 
at At < 200 days in the residuals at wavelengths 
coincident with the 2800A Mgll emission line. 
This decrement is expected based on the results 
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Fig. 14. — Weighted mean of the data-model residu- 
als as a functions of rest wavelength (the bottom-right 
panel of Fig. 13 collapsed onto the wavelength axis). 
The weighted mean of all time lag bins i is defined as 

J2{SFi — SF,nodei,i)wi/J2^i'^^'^^'^ ^ l/f^SFithe Un- 
certainty from Fig. 13. The error bars of the weighted 
means are computed as (^Wi)"^^^. The composite 
quasar spectrum is over-plotted in arbitrary flux {Fx) 
units. The dip at 2800A results from the lagged and 
smeared variability of the Mgll emission line with respect 
to the continuum, and is widened by the errors in esti- 
mating the rest wavelength from the SDSS bands. 



of Reichert et al. (1994), where the Mgll emis- 
sion line is less variable and lags the continuum 
fluctuations by almost 10 days. In Figure 14, 
the time axis is collapsed onto the wavelength 
axis (using the weighted mean) so that the dip 
in the residuals is emphasized. Due to the ex- 
cellent SDSS photometry and large sample size, 
we are able to resolve this feature. Note that the 
large width of the decrement may be due to the 
fact that the rest wavelengths are approximated 
from the fixed effective wavelengths of the SDSS 
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bands. Not accounting for the variation of the ef- 
fective wavelength with the shape and redshift of 
the quasar spectrum sampled limits the accuracy 
of the approximated rest wavelengths. Thus, any 
line feature will be weakened and smoothed by 
the broad bandpasses of the SDSS filters. 

4. Combining Short- and Long-term Quasar 
Variability Measurements (SDSS POSS) 

The exponential distribution analysis in Sec- 
tion 3.2 is sensitive to both SF^o and r through 
the better determined combination defined by 



a = SFoc/^/T. However, for the long-term 
SDSS-POSS data, the Am distribution is mainly 
sensitive to SFoo or r at fixed a. It is then of 
interest to see whether the Am distribution re- 
mains exponential at large \Am\ and if we can re- 
produce the ensemble SF(At) on these long time 
scales. Therefore, we repeat the analysis in Sec- 
tion 3.2 using the 81,189 SDSS quasars that are 
also observed in DPOSS. We only use the DPOSS 
data for the Am analysis because the errors for 
POSS-I are significantly larger. After applying 
the data quality cuts, there are 56,732 total Am 
measurements in the G, R, and / bands. Fig- 
ure 15 shows the distribution of time lags in the 
observed and rest frames. 

The SDSS-DPOSS Am distributions are shown 
in Figure 16. Due to the larger photometric er- 
rors, these distributions appear more Gaussian 
than those restricted to SDSS observations. In 
some panels, the peak is also offset from zero to- 
wards positive Am values. This is likely due to 
the Malmquist-like bias discussed in dV05, where 
many objects that were fainter at the time of the 
DPOSS observation were then lost and not in- 
cluded. In general, the model curves (solid lines) 
are able to reproduce the shape of the distribu- 
tions but underestimate the rms, even when ac- 
counting for the estimated photometric errors of 
CTphot = 0.1 mag. 

Finally, Figure 17 shows the SF over a large 
range of rest-frame time lag, including the SDSS- 
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Fig. 15. — Top: Distribution of epochs for the SDSS and 
DPOSS observations. Middle: Distribution of the SDSS- 
DPOSS time lags in the observer's frame. Bottom: Dis- 
tribution of the SDSS-DPOSS time lags and wavelengths 
in the rest frame. 



POSS data as measured in Ses06. We find an 
overall characteristic time scale of 835 it 82 days 
with SFoo = 0.26 ±0.01 mag. Note that in dV05, 
the SF reaches 0.46 it 0.02 mag at 40 years in the 
rest frame for the g and r bands combined. 
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Fig. 16. — Similar to Fig. 10 but for SDSS-DPOSS observations over the northern DR5 footprint (with i < 19.1 in SDSS). 
The photometric accuracy was assumed to be a phot = 0.1 mag when computing the models. The Malmquist-like bias seen 
in the bottom panels is taken into account by shifting the model Am distribution to the right by the observed mean. 
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Fig. 17. — Top; Similar to Fig. 5, but plotted against rest-frame time lag, for data with rest-frame wavelength in the 
range 2000-3000 A. The data have been corrected for photometric errors by subtracting 0.018\/2 in quadrature. The green 
dotted line shows the prediction of the DRW model trained on S82 data with all r and SFoo values multiplied by 2 and \/2, 
respectively. Bottom: The data are shown along with comparison (;-band fits from VB04 (purple line), Wilhite et al. (2008; 
thin line) and Voevodkin (2011; red dash-triple-dotted line). The normalization for these curves is arbitrary, since we are 
simply comparing the slopes. The blue dash-dotted line shows the best-fit broken power law to our data with the break set 
to 42.3 days, as in Voevodkin (2011). 
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4.1. Data— Model Comparison 

In Figure 17, we show our model SF as de- 
scribed in Section 3.2.1. We again see that the 
model is biased to higher values of the SF be- 
tween 1 and 200 days. We attribute this discrep- 
ancy to one or both of the following: 

1. At short At, the errors in our measured 
SF are underestimated, and some source of 
systematic error, such as an inaccurate cor- 
rection for SDSS photometric errors, causes 
the measured SF to appear smaller than 
what the model predicts. 

2. The DRW model is inaccurate at short 
At and therefore overpredicts the ensem- 
ble SF. However, this is unlikely given that 
the results using Stripe 82 light curves as 
well as densely-sampled OGLE light curves 
(Zu et al. 2012) show the DRW model to be 
a good fit on time scales between ~5 and 
200 days. 

Nevertheless, our measured slope of the SF seems 
well-constrained. We again see that beyond 
200 days, the model is biased to lower values of 
the SF. The difficulty to exactly reproduce the 
long-term variability amplitudes suggests that 
the distributions of long-term amplitudes and /or 
time scales used in our model underestimate the 
true distributions. This bias may result from the 
following: 

1. A bias in the measured r (and thus SFoo) 
values for S82 quasars due to insufficient 
light curve lengths. Since SF{At « t) = 
(Ti/| At|, a = SFoo/\/t is the most strongly 
constrained quantity on time scales shorter 
than r. That is, on a grid of log(r) versus 
log(SFoo), the best-fit values will be scat- 
tered due to fitting errors along lines of con- 
stant (T, but much less perpendicular to it 
(i.e., along lines of constant K = t^/SF^). 
Therefore, as the light curve length de- 
creases, the mean best-fit a will not vary 



significantly, while the best-fit parameters 
r, SFoo, and K will become biased (see 
KozllO; MaclO; MacLeod ct al. 2011). 

2. A bias in the deterministic model (Eq. 6). 
Due to the finite length of the S82 sur- 
vey, 22% of the S82 sample was excluded 
from the analysis in MaclO due to inde- 
terminately long time scales. In this case, 
the measured parameters r and SFoo are 
accurate for the remaining 78% of the S82 
sample, but Eq. 6 is only accurate for the 
lower values of r (r ^ lO"') and SFoo from 
which the correlations are derived. 

In principle, both effects will lead to lower overall 
r and SFoo in the distributions for S82 quasars, 
which will in turn cause the model ensemble SF 
to flatten at AIrf < 700 days as seen here. In- 
deed, the excluded 22% of the sample will mostly 
contribute power to the long-term rather than 
the short-term SF given their indeterminately 
long time scales. However, when simply includ- 
ing an additional population of long-r objects to 
fill in the missing 22%, we are unable to quan- 
titatively reproduce the last SDSS data point at 
AtRF = 1000 days unless all the additional ob- 
jects have r ^ 10, 000 days in the rest frame 
(while keeping the same distribution of a as ob- 
served). On the other hand, if we double all the 
r values in our model (instead of adding an extra 
population of quasars at long r), while keeping 
a fixed so that the SFoo values are multiplied 
by ^/2, we are able to reproduce the observed 
long-term SF, as shown by the agreement be- 
tween the green lines in Figures 17 and 11 to 
the data points at long At^p. In summary, the 
best-fit A coefficients from Eq. 6 need to be al- 
tered upwards by 0.15 dex in the case of SFoo, 
and by 0.30 dex in the case of r, in order to ex- 
plain the long time scale constraints provided by 
the SDSS-POSS dataset. 

At least half of these corrections can be under- 
stood as due to a fitting bias towards shorter r 
(see Figure 7 in MaclO and Figure 15 in MacLeod 
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et al. 2011). Another effect could be due to un- 
certain behavior for long time scales. The DRW 
process corresponds to a power spectral distribu- 
tion (PSD) proportional to 1//^ at frequencies 
/ > (27rr)^^, flattening to a constant at lower 
frequencies. Using the S82 data and computa- 
tional technique described in MaclO, we were 
able to rule out an extrapolation of the 1//^ 
power-law. However, we were unable to distin- 
guish between a 1 j or a 1// PSD at frequencies 
/ < (27rT)~^, where the latter dependence is ob- 
served in the X-ray PSDs for Galactic black holes 
as well as AGN (McHardy et al. 2006; Kelly et 
al. 2011). 

We compare the observed SF slopes from Vo- 
evodkin (2011), VB04, and Wilhite et al. (2008) 
to our data in the bottom panel of Figure 17. 
Voevodkin (2011) found that a broken power-law 
provides a good fit to the S82 (7-band ensemble 
SF with a slope of 0.33, steepening to 0.79 be- 
low 42 days. We also show the best-fit broken 
power law to our data with the break fixed to 42 
days, and we find power law indices at short and 
long At of 0.53 and 0.40, respectively. While our 
2-epoch SDSS data are consistent with the shal- 
lower slope of 0.33, our results do not support the 
conclusion of a much steeper SF(At) for small At 
found by Voevodkin (2011) for either the S82 or 
2-epoch data sets. We interpret the broken power 
law form to be a consequence of the turnover in 
the SF due to the mean characteristic time scale, 
as the observed SF is fully consistent with the 
form expected for a DRW (Eq. 5). 

5. Predictions for Future Surveys 

Using the observed Am distributions, we can 
predict the number of quasars with Am exceed- 
ing an arbitrary limit that might be seen in a 
survey with a given number of quasars. This in- 
formation is useful for transient identification, in 
particular to identify quasars that contaminate 
candidate lists of other objects. For example, 
Vanden Berk et al. (2002) reported an orphan 
gamma-ray burst afterglow based on a 2.5 mag- 



nitude decrease in the optical flux of an uniden- 
tified point source. Instead, as pointed out by 
Gal- Yam et al. (2002), the observations are best 
explained by the presence of a spectroscopically 
identified, highly variable quasar. 

In order to quantify the importance of quasar 
contamination in future transient surveys, such 
as the Palomar Transient Factory (PTF; Law et 
al. 2009) and the Large Synoptic Survey Tele- 
scope (LSST; Ivezic et al. 2008), we need to know 
the probability that a quasar can increase bright- 
ness by Am (over the faint survey limit) within a 
time At. For example, assume that no source is 
detected above a faint limit of m faint on a given 
night, but when repeating the observation some 
time later, a source is detected with magnitude 
m = m faint ~ Am. In this case we would like to 
know how many quasars with a Am at least as 
large as observed could be present in a particu- 
lar scanned area and to a faint limit of m faint ■ 
Given that our stochastic model performs well at 
reproducing the observed ensemble variability of 
quasars, we can make robust and useful predic- 
tions for future surveys. 

We make use of a mock LSST quasar sam- 
ple which includes absolute B magnitudes {Mb) 
and redshifts generated over 100 deg2 of sky us- 
ing the luminosity function from Bongiorno et al. 
(2007) for the purposes of LSST image simula- 
tions. We limit the sample to Mb < — 20 follow- 
ing Table 10.2 in the LSST Science Book (LSST 
Science Collaborations and LSST Project 2009). 
The distance modulus is computed assuming a 
standard cosmology: = 0.27, JIa = 0.73, 

and h = 0.71. The LSST magnitudes (m) and 
rest-frame Mj (based on Mb) are computed us- 
ing the composite quasar spectrum from Vanden 
Berk et al. (2001). The black hole masses are 
estimated from the Mj values using the prescrip- 
tion in MaclO. 

We generate magnitude differences for each 
quasar in the mock LSST sample based on the 
DRW model most appropriate for the quasar's 
physical parameters, as described in Section 3.2.1. 
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When computing the DRW model, the r (SFqo) 
values are multiplied by 2 {V^) with respect to 
the expected values based on MaclO in order to 
correct for the bias due to limited time sampling 
in S82 (see previous section). We note that cor- 
recting for this bias should also account for at 
least some of the residual scatter in K estimated 
in Section 3.2.1; however, for simplicity we still 
include this scatter along with the bias correc- 
tion in our simulations. First, we consider three 
different survey faint limits of mjaint = 19.1 
(to establish similarity with the SDSS results), 
m faint = 22, and m faint = 24.5, excluding ah 
mock quasars with m > m faint in each case. For 
the m < 19.1, 22, and 24.5 simulations respec- 
tively, 1000, 100 and 11 model Am values are 
generated per quasar to increase the sample size. 
A Gaussian noise component of width V^aphot 
is added to all curves to simulate a photometric 
accuracy similar to the SDSS (the photometric 
accuracy for future surveys such as LSST will 
likely be better than that for the SDSS, but this 
is a higher order question than investigated here) . 
Note that we retain all Am values in the simula- 
tions, including \Am\ > 3. 

Figure 18 shows the simulated cumulative dis- 
tribution of Am in the urz bands for three faint 
magnitude limits. Also shown are Gaussian ana- 
lytic functions as thin curves with the same rms 
as the data. Using these cumulative distribu- 
tions, which are based on ~10^ mock quasars, 
we predict probabilities down to our resolution 
limit (P ;^ lO"''). Table 1 lists the predicted 
probabilities of observing a quasar with m < 24.5 
and a magnitude difference of Am > 1 mag and 
Am > 2 mag over 3, 30, and 300 days, us- 
ing the observed frame, in the urz bands. Due 
to the exponential nature of the Am distribu- 
tions, the probability of observing Am > 1 mag 
reaches 0.02 in the u band (where variability is 
strongest) for time lags of 300 days, and 6 x 10~^ 
for Am > 2 mag. Assuming Gaussian Am dis- 
tributions will result in erroneous probability es- 
timates of 9 X 10~^ and 10~^, respectively. 
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Fig. 18. — The predicted cumulative distribution of 
magnitude difTerences in u, r, and z bands as a func- 
tion of observed time lag (At) for ~1 million quasars with 
magnitudes less than 19.1 (top three panels), 22 (middle 
three panels), and 24.5 (bottom three panels). The thin 
curves show the Gaussian distributions with the same rms 
widths. 
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Fig. 19. — The predicted cumulative distribution of g- 
band (top) and r-band (bottom) magnitude differences, 
including realistic PTF photometric errors as a function 
of observed time lag (At) for quasars brighter than the 
PTF magnitude limit [g < 21.0, r < 20.6). The thin 
curves show the Gaussian distributions with the same rms 
widths. 

Next, we adopt realistic PTF photometric er- 
rors which vary as a function of magnitude, and 
repeat the simulation in g and r using the PTF 
magnitude limits (r < 20.6, g < 21.0). The er- 
rors range from 0.008 mag at the brightest lim- 
its to 0.2 mag at the faint end. As shown in 
Figure 19, the resulting distributions are more 
Gaussian than those in the previous simulations 
due to the different dependence of errors on mag- 
nitude. However, at extreme values of Am, the 
Gaussian curves significantly underestimate the 
detection probabilities in both g and r. 

These results would have been useful to Van- 
den Berk et al. (2002), who reported a single 
transient with Am = 2.5 mag over 410 days 
in the gri bands in early SDSS data covering 
1, 500 deg^. Since the object had the spectrum of 
a normal galaxy in its faint phase (rather than a 
nonstellar spectrum with broad emission lines in- 
dicative of quasars), the authors concluded that 
it was more likely to be a gamma-ray burst af- 
terglow than a highly variable quasar, given the 
large drop in flux. We find that the probability 



of a quasar having Am = 2.5 mag over 410 days 
with m faint = 22 is 10~^ in the g and r bands. 
To find the number of quasars expected to ex- 
hibit this variability in 1,500 deg^ of sky, one 
needs to know the quasar density. We extrap- 
olate Figure 13 in Richards et al. (2006) to a 
density of ~120 deg~^ at i < 22. Therefore, 
we expect roughly 10^ quasars in 1,500 deg^ to 
i < 22. Given a probability of 10^^, we expect 
that roughly one quasar will show Ar = 2.5 mag 
over At = 410 days in the SDSS sample studied 
by Vanden Berk et al. (2002). This is consistent 
with their finding, and with the follow-up study 
of Gal- Yam et al. (2002) who confirmed that it is 
indeed a highly variable quasar. If one assumes 
a Gaussian distribution of Am, the probability 
is less than 10~^, and one would expect ^0.1 
quasar to be found with these parameters. Given 
this expectation value, the probability to detect 
one quasar is at most 0.09, and the quasar hy- 
pothesis can be (erroneously) rejected at a ;^ 2o" 
level. 

6. Discussion 

We have assembled, organized and publicly 
released a dataset including ~3.5 million pho- 
tometric measurements for 80,000 spectroscopi- 
cally confirmed quasars. The available time lags 
span 0.8 days to almost 20 years in the observer's 
frame. We have analyzed and quantified the 
observed variability in the observer's and rest 
frames. By assuming a DRW model for each 
quasar in our sample, we reconcile the observed 
variability of individual quasars in S82 with their 
ensemble statistics. Our principal results are as 
follows: 

1. Long-term quasar variability measurements, 
constrained using SDSS and POSS data 
for time lags up to 50 years (in the ob- 
server's frame), conclusively show that a 
simple power-law dependence for the struc- 
ture function cannot be extrapolated be- 
yond a decade, and suggests an average 
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Table 1 

Predicted Am Probabilities for m < 24.5 



At (days) 


Band 


P(Am > 1 mag) 


P(Am > 2 mag) 
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< 2 X 10-6 


< 2 X 10-6 
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< 2 X 10-6 


< 2 X 10-6 
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< 2 X 10-6 


< 2 X 10-6 


30 
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7 X 10-^ 


< 2 X 10-6 


30 
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4 X 10-6 


< 2 X 10-6 


30 
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< 2 X 10-6 


< 2 X 10-6 


300 


u 


0.02 


6 X 10-^ 


300 


r 


0.005 


6 X 10-^ 


300 


z 


0.002 


1 X 10-5 



characteristic time scale for quasar variabil- 
ity in the rest frame of ~2 years and an 
average long-term dispersion of ^0.26 mag 
(for rest wavelengths 2000-3000A). This 
behavior extrapolates well to the UV re- 
sults of Welsh, Wheatley & Neil (2011), 
who find that the SF for GALEX NUV 
data reaches about 0.4 mag and flattens at 
AtRF > 300 days. This SF hmit corre- 
sponds to a limiting SF value of 0.33 mag 
when using our definition of the SF (see 
Section 3.1.2), and 0.27 mag when also scal- 
ing to the «-band using the wavelength de- 
pendence from Eq. 6. This result is in 
close agreement with the SF at the short- 
est wavelengths in our data set. Voevodkin 
(2011) found that a broken power-law pro- 
vides a good fit to the S82 ensemble SF 
with a slope of 0.33 at long time scales 
steepening to 0.79 below 42 days. Our 2- 
epoch SDSS data are consistent with the 
shallower slope of 0.33, but our data do not 
support the conclusion of a much steeper 
SF(At) for small At found by Voevodkin 
(2011). While we cannot rule out a broken 
power-law dependence with the available 



data, the observed SF is fully consistent 
with the form expected for a DRW (Eq. 5). 

2. We tested the DRW model results based 
on SDSS Stripe 82 data on an indepen- 
dent dataset, and confirm that the vari- 
ability parameters r and SF^o correlate 
with physical parameters as found for indi- 
vidual quasars (e.g., MaclO and references 
therein). This is evident from the agree- 
ment of our model with the observed en- 
semble variability of SDSS quasars. How- 
ever, the results indicate that the measured 
T and SFqo distributions are biased low 
for the S82 sample by a factor of about 
2 and \/2) respectively. This bias most 
likely results from the 10-year length limit 
of the S82 light curves, although it could 
also be due to uncertain behavior for long 
time scales. The best-fit A coefficients from 
Eq. 6 (MaclO) need to be shifted upwards 
by 0.15 dex in the case of SFoo, and by 
0.30 dex in the case of r, in order to ex- 
plain the long time scale constraints pro- 
vided by the SDSS-POSS dataset. These 
shifts leave the shorter time scale variabil- 
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ity statistics unchanged. 

3. For a given time lag and wavelength, the 
magnitude difference (Am) distribution is 
exponential rather than Gaussian for large 
magnitude changes. This is well explained 
as a cumulative effect of averaging over 
quasars with a range of different r and 
SFoo- This is a remarkable result given 
that the Am distribution of every individ- 
ual quasar is Gaussian. 

4. We made predictions for the incidence of 
quasar contamination in transient surveys 
using detailed simulations of quasar light 
curves from a mock LSST catalog. Due to 
the exponential nature of the Am distri- 
butions for quasars, the probability of ob- 
serving Am > 1 mag reaches 0.02 in the 
u band (where variability is strongest) for 
time lags of 300 days, and 6 x 10""^ for 
Am > 2 mag. Assuming Gaussian Am 
distributions will result in erroneous likeli- 
hood estimates that are about 10 and 1000 
times smaller, respectively. 

It is clear that a major limitation for the S82 
quasars is the quality of light curves in both sam- 
pling density and time span. It is also clear that 
our variability model needs to be better tested 
given the evidence for a likely bias in the S82 
time scale estimates. The best current sample 
for these improvements is that from the OGLE 
microlcnsing survey, since the light curves are 
more densely-sampled and longer than for S82. 
Here, the probem is the lack of spectroscopic 
identification of quasar candidates, although the 
follow-up confirmation of quasars is rapidly im- 
proving (Kozlowski et al. 2011b). The next- 
generation surveys will also greatly improve the 
constraints on the long-term SF both individ- 
ually and for ensembles of quasars. The best 
short-term prospects are PTF (Law et al. 2009), 
Pan-STARRS (Kaiser et al. 2002), and the Dark 
Energy Survey (DES; Honscheid et al. 2008). 
In particular, the DES supernova program will 



greatly expand many of the S82 quasar light 
curves with griz sampling once per week for 
~3 months per year over 5 years. The combi- 
nation of SDSS, Pan-STARRS, DES, and LSST 
will yield well-sampled light curves covering over 
25 years for 10,000 quasars in Stripe 82. 

The success of the model presented here sug- 
gests that a range of characteristic time scales 
exists among an ensemble of quasars, which can 
be related to physical time scales in the accretion 
disk. While we assumed a single r per quasar, 
there is evidence that multiple time scales can 
exist for a given quasar (Collier & Peterson 2001; 
Kelly et al. 2011). Therefore, the study presented 
here can be extended to adopt the model in Kelly 
et al. (2011), which fits more than one r for a 
given object. This model, also called a mixed 
Ornstein-Uhlenbeck (OU) process, reproduces a 
PSD of the form exhibited by the X-ray light 
curves of galactic black holes and AGN, which 
is flat below a low-frequency break, decays as 
1/f above the low- frequency break, and steep- 
ens to 1//^ above a high-frequency break. In 
this case, with two characteristic time scales for 
each quasar, the long-term ensemble SF can be 
revisited and possibly explained in the context 
of a mixed OU process. Note that recent op- 
tical data from the Kepler mission (Mushotzky 
et al. 2011), which have a sampling of 1 data 
point roughly every 30 minutes and 0.1% errors, 
suggest an additional break to a steeper slope 
(~ 1//^), but this dependence is seen on time 
scales shorter than can be resolved in SDSS data. 
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Appendix 

Here we present a database of ~3.5 million 
photometric measurements for 80,000 spectro- 
scopically confirmed quasars to be used for time 
variability studies. The database consists of 
three different data sets: two with repeated SDSS 
imaging in five UV-to-IR photometric bands, and 
one with SDSS versus POSS imaging for three 
bands. The observed time lags span the range 
from 0.8 days to 10 years for the SDSS data sets, 
and up to 20 years for the SDSS vs. DPOSS data 
set. The three data sets are described below: 

1. Northern Survey: The SDSS imaging data 
are obtained by drift scanning. Because 
of the scan overlaps, and because of the 
scan convergence near the survey poles, 
about 40% of the northern survey area 
(~4000 deg^) is surveyed at least twice. 
This method provides 2-epoch 5-band cov- 
erage for ~25,000 spectroscopically con- 
firmed quasars. We adopt the SDSS BEST 
photometry listed in the DR7 Quasar Cat- 
alog V (Schneider et al. 2010) for the pri- 
mary observations, and we searched for all 
unresolved secondary observations within 1 
arcsecond of the primary using CasJobs. 
Hero, wo include all observations up to 



This 2-column preprint was prepared with the AAS IM^iX 
macros v4.0. 
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three per quasar (only 0.2% of the sam- 
ple had more than three observations). We 
also include the redshifts, absolute magni- 
tudes, FIRST, RASS, and 2MASS photom- 
etry as listed in the DR7 Quasar Catalog 
V, and the black hole masses as measured 
from emission line widths by Shen et al. 
(2011). The catalog format is found in Ta- 
ble 2. We also provide a list of all the ob- 
jects (in the same order) with the same ex- 
act format as the DR7 Quasar Catalog V 
(Table 2 in Schneider et al. 2010). 

2. Southern Survey: About 290 square de- 
grees of the southern survey area has al- 
ready been observed ~60 times to search 
for variable objects and, by stacking the 
frames, to go deeper. This is the SDSS 
Stripe 82, which is 22h 24m < R.A. < 04h 
08m and |Dcc| < 1.27 deg. These multi- 
epoch data have time scales ranging from 
3 hours to almost 10 years. This method 
provides well-sampled 5-band light curves 
for an unprecedented number of quasars 
(9,258). The catalog format is found in 
Table 3, and the light curve file format is 
found in Table 4. We also provide a list of 
all the objects (in the same order) with the 
same exact format as the DR5 Quasar Cat- 
alog IV (Table 2 in Schneider et al. 2007). 

3. SDSS-DPOSS: We also include a catalog 
of all SDSS DR7 quasars with DPOSS ob- 
servations. Following the procedure out- 
lined in Sesar et al. (2006), we have re- 
calibrated DPOSS data (Djorgovski et al. 
1998) in 8,000 deg^ of sky from the SDSS 
Data Release 5. The main advantage of this 
data set, which includes 81,189 quasars, is 
its long time baseline of 20 years. Here, 
we present the SDSS-DPOSS photometry 
in GRI bands, where the latter is accurate 
to 0.10-0.15 mag. The catalog format is 
presented in Table 5. 



For more details, see 

http://www.astro.washington.edu/users/ivezic/macleod/qso_dr7/. 



29 



Table 2 

Northern Survey: Catalog Format 



Column 


Format 


Units 


Description 


1 


16 




Row of DR7 Quasar Catalog V (out or 105783) 


2 


F7.3 


mag 


Absolute i band magnitude'' 


3 


F6.4 




Redshift 


4 


F6.3 


[Mq] 


Black hole mass ("0" means none could be found) * 


5 


F6.3 


[ergs/s] 


Bolometric luminosity^ 


6 


F6.3 


mag 


Galactic extinction m u niter (from ocnlegel et al. 1998) 
b IRb i Peak 20 cm flux density 


7 


F7.3 


mag 


8 


F8.3 




bignal-to- noise ratio tor t IRS i flux density 


9 


F8.3 


[ct/sj 


RAbb ribO/rbO full band count rate 


10 


F7.3 




bignal-to-noise ratio tor RAbb count rate 


11 


F7.3 


mag 


2MASS J band magnitude 


12 


F6.3 


mag 


Error in J 


13 


F7.3 


mag 


2JVlAbb ii band magnitude 


14 


F6.3 


mag 


Error in H 


15 


F7.3 


mag 


2MASS K band magnitude' 


16 


F6.3 


mag 


Error in K 


17 


13 




Total number of observations 


18 


11 




Flag for Stripe 82 object'' 


19 


F9.3 


days 


Modified Julian Date of imaging observation 


20 


F6.3 


mag 


BEST SDSS u band PSF magnitude'' 


21 


F6.3 


mag 


Error in u 


22 


F6.3 


mag 


BEST SDSS g band PSF magnitude'' 


23 


F6.3 


mag 


Error in g 


24 


F6.3 


mag 


BEST SDSS r band PSF magnitude'' 


25 


F6.3 


mag 


Error in r 


26 


F6.3 


mag 


BEST SDSS i band PSF magnitude'' 


27 


F6.3 


mag 


Error in i 


28 


F6.3 


mag 


BEST SDSS z band PSF magnitude'' 


29 


F6.3 


mag 


Error in z 

Fields 19 through 29 repeated for each additional observation 



*A value of 1 (0) indicates a (non-)Stripe 82 source. 

''SDSS photometric measurements arc asinh magnitudes (Lupton, Gunn, & Szalay 1999) and are normalized (to ~3% accuracy) to 
the AB-magnltude system (Oke & Gunn 1983). Uncorrected for Galactic extinction. A value of 0.000 Indicates that the value could 
not be retrieved from the SDSS database. 

''Where Ag, Ar, Ai, A^ = 0.736, 0.534, 0.405, 0.287 x Au, respectively 
''in AB magnitudes; -2.5 log(/„/3631Jy). 

''X-ray data from ROSAT All Sky Survey Bright and Faint source catalogs. A value of —9.000 Indicates a non-detection. 

^All 2MASS data are from the 2MASS All-Sky Data Release Point Source Catalog (PSC) as of 2003 March 25. Note that 2MASS 
measurements are Vega-based, not AB, magnitudes. A value of 0.000 indicates a non-detection. 

^Taken from the Shen et al. (2011) catalog. 

''For Clm = 0.300, CIa = 0.700, h = 0.70, and ag = -0.50. 
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Table 3 

Southern Survey: Catalog Format 



Column 


Format 


Units 


Description 


1 


17 




The name of the hght curve file 


2 


F10.6 


deg 


Median Right Ascension in decimal degrees (J2000) 


3 


F10.6 


deg 


Median Declination in decimal degrees (J2000) 


4 


15 




Row of DR5 Quasar Catalog IV (out of 77429) 


5 


F7.3 


mag 


Absolute i-band magnitude, K-corrected to z = 0'^ 


6 


F7.3 


mag 


Absolute j-band magnitude, K-corrected to « = 2" 


7 


F6.4 




Redshift 


8 


F5.3 


[Mo] 


Black hole mass ("0.000" means none could be found)'' 


9 


F6.3 


[ergs/s] 


Bolometric luminosity'' 


10 


F6.3 


mag 


SDSS BEST u-band PSF magnitude'' 


11 


F6.3 


mag 


SDSS BEST g-band PSF magnitude'' 


12 


F6.3 


mag 


SDSS BEST r-band PSF magnitude'' 


13 


F6.3 


mag 


SDSS BEST t-band PSF magnitude'' 


14 


F6.3 


mag 


SDSS BEST z-band PSF magnitude'' 


15 


F6.3 


mag 


Galactic extinction in u filter (from Schlegel et al. 1998)"^ 



^Taken from the Shen et aL (2008) catalog. A value of indicates a newly confirmed DR7 quasar, and thus is not present in Shen 
et al. (2008) (sec the latest version, Shen et al. 2011, for the values). 

''SDSS photometric measurements are asinh magnitudes (Lupton, Gunn, & Szalay 1999) and are normalized (to ~3% accuracy) to 
the AB-magnitude system (Okc & Gurm 1983). Uncorrected for Galactic extinction. A value of 0.000 indicates that the value could 

not be retrieved from tlie SDSS database. 

'^Where Ag, Ar, Ai, A^ = 0.736, 0.534, 0.405, 0.287 x A„, respectively. This value is set to zero if it is a newly confirmed DR7 quasar 
(see the DR7 Quasar Catalog V for the true values). 

''For = 0.300, Ha = 0.700, h = 0.70, and aq = -0.50. 
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Table 4 

Southern Survey: Light Curve File Format 



Column 


Format 


Units 


Description 


1 


D 


days 


Modified Julian Date for u-band observation 


2 


F6.3 


mag 


Apparent w-band magnitude'' 


3 


F5.3 


mag 


Error in u 


4 


D 


days 


Modified Julian Date for g-band observation 


5 


F6.3 


mag 


Apparent gr-band magnitude" 


6 


F5.3 


mag 


Error in g 


7 


D 


days 


Modified Julian Date for r-band observation 


8 


F6.3 


mag 


Apparent r-band magnitude" 


9 


F5.3 


mag 


Error in r 


10 


D 


days 


Modified Julian Date for i-band observation 


11 


F6.3 


mag 


Apparent t-band magnitude" 


12 


F5.3 


mag 


Error in i 


13 


D 


days 


Modified Julian Date for z-hand observation 


14 


F6.3 


mag 


Apparent z-band magnitude" 


15 


F5.3 


mag 


Error in z 


16 


D 


deg 


Median Right Ascension in decimal degrees (J2000)'' 


17 


D 


deg 


Median Declination in decimal degrees (J2000) 



"Not corrected for Galactic absorption. Bad observations are printed as "-99.99". 
''360 deg is subtracted from all RA values exceeding 300 deg. 
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Table 5 

SDSS-DPOSS: Catalog Format 



Column 


Format 


Units 


Description 


1 


FIG. 6 


deg 


"r-\ • 1 1 A • 'I'll / ~f c\r\r\r\ \ 

KigJit Ascension m decimal degrees (J 2000) 


2 


F10.6 


deg 


"r-\ !• J • "I'll / ic\r\r\r\\ 

Decimation m decimal degrees (J 2000) 


3 


13 




DPOSb Plate ID 


4 


r 9.4 


yrs 


iivpocn oi trie UPUoo observation m tne Cr bana 


5 


F9.4 


yrs 


ijjpGcn or the DrUbb observation in the ti band 


6 


F9.4 


yrs 


Epoch or the DFObo observation m the I band 


7 


F6.2 


mag 


DFOsb catalog G magnitude 


8 


F6.2 


mag 


Error in G 


9 


F6.2 


mag 


G-band quality nag 


10 


F6.2 


mag 


DPOSS catalog it magnitude 


11 


F6.2 


mag 


Error in ii 


12 


F6.2 


mag 


it- band quality nag 


13 


F6.2 


mag 


DFUbo catalog i magnitude 


14 


F6.2 


mag 


Error in J 


15 


11 


mag 


7-band quality flag° 


16 


16 




Row of DR7 Quasar Catalog V (out of 105783) 


17 


F7.3 


mag 


Absolute i band magnitude*" 


18 


F6.4 




rcdshift 


19 


F6.3 


[Mo] 


Black hole mass ("0.000" means none could be found)^ 


20 


F6.3 


[ergs/s] 


Bolornctric luminosity^ 


21 


F6.3 


mag 


Galactic extinction in u filter (from Schlegel et al. 1998) 


22 


11 




Number of observations 


23 


11 




Flag for Stripe 82 object'' 


24 


15 


days 


Modified Julian Date for SDSS imaging observation 


25 


F6.3 


mag 


BEST SDSS g band PSF magnitude'' 


26 


F6.3 


mag 


Error in g 


27 


F6.3 


mag 


BEST SDSS r band PSF magnitude'' 


28 


F6.3 


mag 


Error in r 


29 


F6.3 


mag 


BEST SDSS i band PSF magnitude'' 


30 


F6.3 


mag 


Error in i 



^A value of 1 (0) indicates a (non-)Stripe 82 source. 

''SDSS photometric measurements are asinh magnitudes (Lupton, Gunn, &; Szalay 1999) and are normalized (to ~3% accuracy) to 
the AB-magnitude system (Oke & Gunn 1983). Uncorrected for Galactic extinction. A value of 0.000 indicates that the value could 

not be retrieved from the SDSS database. 

''Where Ag, Ar, A,, A^ = 0.736, 0.534, 0.405, 0.287 x Au, respectively. 

''Not corrected for Galactic absorption. Bad observations are printed as "-99.99". 

Tlag indicating whether the recalibrated magnitude is of good quality (a value of 1 is "good"). See Ses06 for a description of these 
cuts. 

^Taken from the Shcn ct al. (2011) catalog. 

■Tor Q.m = 0.300, Oa = 0.700, h = 0.70, and ag = -0.50. 
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